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Abstract. The global stratospheric ozone-layer depletion results in an increase in 
biologically harmful ultraviolet (UV) radiation reaching the surface and penetrating to 
ecologically significant depths in natural waters. Such an increase can be estimated on a 
global scale by combining satellite estimates of UV irradiance at the ocean surface from 
the Total Ozone Mapping Spectrometer (TOMS) satellite instrument with the SeaWiFS 
satellite ocean-color measurements in the visible spectral region. In this paper we propose 
a model of seawater optical properties in the UV spectral region based on the Case 1 
water model in the visible range. The inputs to the model are standard monthly SeaWiFS 
products: chlorophyll concentration and the diffuse attenuation coefficient at 490nm. 
Penetration of solar UV radiation to different depths in open ocean waters is calculated 
using the RT (radiative transfer) quasi-single scattering approximation (QSSA). The 
accuracy of the QSSA approximation in the water is tested using more accurate codes. 
The sensitivity study of the underwater UV irradiance to atmospheric and oceanic optical 
properties have shown that the main environmental parameters controlling the absolute 
levels of the UVB (280-320nm) and DNA-weighted irradiance underwater are: solar- 
zenith angle, cloud transmittance, water optical properties, and total ozone. Weekly maps 
of underwater UV irradiance and DNA-weighted exposure are calculated using monthly- 
mean SeaWiFS chlorophyll and diffuse attenuation coefficient products, daily SeaWiFS 
cloud fraction data, and the TOMS-derived surface UV irradiance daily maps. The final 
products include global maps of weekly-average UVB irradiance and DNA-weighted 
daily exposures at 3m and 10m, and depths where the UVB irradiance and DNA- 
weighted dose rate at local noon are equal to 10% of their surface values. 

1. Introduction 

Increased levels of biologically harmful UV-B radiation (280-320nm) resulting from the 
depletion of Earth’s ozone layer [WMO, 1999; Herman et al., 1996] have been shown to 
affect aquatic ecosystems. One of the important effects of enhanced levels of UVB 
radiation is a reduction in the productivity of phytoplankton caused by their moving to 
increased depths where there is less PAR (photosynthetic active radiation). This effect 
has been conclusively confirmed in the laboratory [ Worrest and Hader, 1989] and from 
field observations in the Bellinghausen Sea (Antarctica). Smith et al. [1992] compared the 
productivity of Antarctic phytoplankton beneath and outside the Antarctic ozone hole and 
reported a 6-12% loss due to enhanced UVB. Other studies report smaller effects, 
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generally ranging from 1% [Arrigo, 1994] to 5% [Neale et al. 1998a, b; Holm-Hansen et 
al. 1993]. Enhanced UVB radiation could also affect the photochemical production of 
carbonyl sulfide (COS) in seawater [ Zepp and Andreae, 1994], thereby augmenting the 
greenhouse effect and affecting other long-term global biogeochemical cycles. 

The quantitative assessment of UV effects on aquatic organisms on a global scale 
requires an estimate of the in-water radiation field. The total ozone and UV reflectivity 
measurements from the Total Ozone Mapping Spectrometer (TOMS) satellite 
instruments allow calculation of global daily UV irradiance at the ocean surface [Eck et 
al. 1995; Herman et al. 1996; Krotkov et al. 1998; Herman et al. 1999]. Because of long 
time record (since 1978) and contiguous global spatial coverage, TOMS data are vital not 
only for estimating zonal average trends in surface UV irradiance [Madronich 1992; 
Herman et al. 1996], but also in studying geographical differences in surface UV 
climatology [Herman et al., 1999; WMO 1999]. The current TOMS instrument onboard 
NASA’s Earth-Probe satellite has been collecting data since July 1996, and is expected to 
continue ozone and surface UV reflectivity measurements to provide overlap with launch 
of the last TOMS instrument (July 2000). Additional spacebome ozone sensors are 
planned by the United States and other countries well through the next century (Global 
Ozone Monitoring Experiment (GOME) on European ERS-2 satellite, Triana-EPIC 
(Earth Polychromatic Imaging Camera), and Ozone Monitoring Instrument (OMI) on 
NASA’s EOS-CHEM satellite). 

Estimates of UV transmission in ocean surface waters require knowledge of the inherent 
and apparent optical properties of seawater. For ocean properties, the Coastal Zone Color 
Scanner flown onboard NASA's Nimbus-7 satellite and current ocean-color satellite 
instruments, such as NASA's SeaWiFS instrument (Sea-viewing Wide Field-of-view 
Sensor, Hooker et al., 1992) were designed to provide frequent global measurement of 
the water-leaving radiances from which inherent water optical properties (IOP) and 
constituents (Chlorophyll concentration) can be inferred. 

The main goal of this paper is to estimate the penetration of UV light into Case I open 
ocean waters using large-scale TOMS surface-UV and SeaWiFS ocean-color 
measurements. In assimilating these satellite data sets, two major problems arise: the fast 
radiative transfer modeling of the penetration of UV light into the water and extrapolation 
of water optical properties derived from the satellite visible channels to the UV spectral 
region. The paper discusses both problems with focus on the "radiative transfer" (RT) 
part. In section 2 we briefly discuss the SeaWiFS and TOMS data as input to the RT 
models. Section 3 discusses the water RT model in more details as well as the 
parameterization of the UV optical properties for Case I water. The accuracy of the water 
RT approximation is tested using well-established accurate models ( Mobley et al., 1993). 
Section 4 discusses the calculated maps of underwater UV irradiance including the depth 
dependence of the UV dose rates for UVB and unshielded DNA. The sensitivity study to 
water optical properties is also discussed. 
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2. Satellite Data 


The general algorithm (data flow) for incorporating both TOMS and SeaWiFS satellite 
data into radiative transfer models is shown in Figure 1. In section 2.1 we discuss the 
SeaWiFS data in more detail. Section 2.2 discusses TOMS surface-UV product. 

2.1 SeaWiFS data 

We use the Level 3 binned monthly mean SeaWifs data for estimates of chlorophyll concentration and the 
diffuse attenuation coefficient A^(490nm), as well as the daily SeaWiFS cloud fraction data. The 
SeaWiFS instrument has 8 bands between 412 nm and 865 nm [Barnes, et al., 1994], and 
is in descending (North to South) noontime sun-synchronous orbit. The orbit and swath 
characteristics allow global coverage every two days, with gaps between consecutive 
swaths at low latitudes being filled the following day. The inherent resolution of the 
instrument is about 1 km at nadir, but the coverage is subsampled on-board the spacecraft 
to provide 4 km global area coverage (GAC). The global coverage includes both land 
and ocean areas, but is limited to solar zenith angles less than 70° to conserve spacecraft 
power. The GAC data is processed at the SeaWiFS Project Office in near-realtime. 

The stability of the sensor is tracked using a time series of monthly lunar calibration 
images [Barnes, et al., 1999]. The prelaunch laboratory calibrations [Johnson, et al., 
1999] are adjusted so that the derived water-leaving radiances in bands 1-6 (412 nm - 670 
nm) are matched to field measurements collected by the Marine Optical Buoy (MOBY) 
located off Lanai, Hawaii [McClain, et al., 1998]. The calibrated radiances (Level-1) 
over the ocean are atmospherically corrected [Gordon and Wang, 1994] to derive Level-2 
geophysical products, e.g., normalized water-leaving radiances [Gordon and Clark, 
1982], chlorophyll-a [O’Reilly, et al., 1997], and diffuse attenuation coefficient at 490 
nm, ATrf(490), [Mueller and Trees, 1997]. These data are spatially binned and averaged on 
a 9 km global grid (Level-3) for each day. The daily gridded data is temporally averaged 
at 8-day, monthly, and annual periods. Figures 2 and 3 show the July 1998 monthly- 
average chlorophyll and ^(490) fields. McClain, et al. [1998] presented the initial 
validation of the chlorophyll retrievals in Case 1 waters, i.e., waters where the reflectance 
is determined by absorption. Further analyses using a larger data set supports the initial 
evaluations that the chlorophyll retrievals are within the SeaWiFS accuracy goal (± 35% 
for Case 1 waters) [Hooker and McClain, 1999]. 

For underwater irradiance calculations, one needs to know both the direct and diffuse 
components of the surface irradiance, and the boundary conditions at the air-water 
interface. For ocean water, the boundary conditions are largely determined by wind 
speed. The TOMS standard data (described below) provides only the global surface 
irradiance (diffuse plus direct). To calculate the daily-average direct irradiance, 
information on average cloud fraction in each grid-cell is required. Such information is 
obtained from the 865nm channel of the SeaWiFS sensor. 

The SeaWiFS cloud-albedo threshold over ocean was set at 1.1% [Arrigo and McClain 
1995]. The relation is a binary one: if the threshold is crossed, the SeaWiFS GAC pixel 
(4km by 4km) is declared cloud contaminated and the cloud flag is set for that pixel. 


3 



During binning, if the flag is set, that sample is considered to be 100% cloud, otherwise 
0%. The threshold albedo of 1.1% agrees with model calculations of reflected radiance 
for spectral band 0.82-0. 88pm above the ocean for a clear atmosphere (the normalized 
reflection coefficient (albedo) is in the range 0.5-1% for clear Rayleigh atmosphere and 
2- 4% including haze and thin cirrus clouds effects). The current cloud flag also masks 
sun glint, high aerosols, and thin cirrus clouds. 

Figure 4 shows a 5-day average SeaWiFS cloud fraction map averaged over 1 by 1 deg 
grid cells. Because of the conservative cloud flag threshold required by the SeaWiFS 
ocean color algorithm, there are relatively small amounts of clear regions (Mediterranian, 
Tropical Atlantic Region, and Equatorial Pacific belt). On the other hand, regions 
declared to have 100% cloud cover may actually contain aerosol plumes that originated 
on land (Sahara Region) and were transported over open ocean regions [Herman et al., 
1997], 

2.2 TOMS data 

TOMS measures solar radiation backscattered from the Erath’s atmosphere in 6 narrow 
(~lnm) wavelength channels (308, 312, 317, 322, 331, and 360nm for Earth-Probe, July 
1996 to present, and 312, 317, 331, 340, 360, and 380 nm for Nimbus-7, November 1978 
to April 1993). This information is converted into total column ozone amounts, scene 
reflectivity, aerosol amounts, and cloud transmittance [McPeters et al., 1996]. The Level 
2 TOMS products are obtained for areas corresponding to the satellite field of view 
(FOV) on the Earth’s surface (~40 km by 40km in nadir direction for EP/TOMS since 
January 1998 to the present). Since the TOMS orbit is Sun synchronous, the FOV can be 
combined into a global map of the same parameters at the same solar time every day 
(Level 3 products). Because of the accurate in-flight calibrations and long-term stability 
of the instrument, the TOMS satellite products are able to maintain absolute and relative 
accuracy over wide-geographical areas and long periods. 

The TOMS daily gridded (Level 3) products (ozone, reflectivity and aerosol index) are 
used as an input for a radiative transfer model to generate daily global maps of the 
surface global (direct plus diffuse) spectral irradiance at the satellite overpass time 
[Krotkov et al., 1998; Herman et al., 1999]. To calculate daily UV exposures we use 
cloud and aerosol amounts at the overpass time, but neglect their diurnal variation. 
Detailed information about TOMS products and daily erythemal maps of the surface LTV 
can be found on the TOMS web site: http://toms.gsfc.nasa. gov and in the TOMS User’s 
Guide [McPeters et al., 1998]. 

The accuracy of estimated UV exposure (time integrated UV irradiance) at a specific 
location is limited by the coarse spatial resolution of the TOMS observations and the 
satellite’s ability to obtain only one UV irradiance estimate per day (except for high 
latitudes). Because of the highly variable nature (temporal and spatial) of cloud and 
aerosol cover the TOMS daily UV estimations should be averaged over periods of at least 
a week to obtain a good estimate of the accumulated UV exposure [Herman et al., 1999]. 
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For the purposes of estimating UV irradiance at the Earth’s surface, there are two classes 
of aerosols that must be considered: 1) aerosols that only scatter UV radiation and 2) 
aerosols that both scatter and absorb UV radiation. The first category is included in the 
TOMS measured scene reflectivity, and attenuates UV radiation in a manner that very 
closely approximates clouds of equivalent reflectivity. Attenuation values of up to 10% 
are common. 

The second category, absorbing aerosols, must be treated separately. First, clouds are 
distinguished from aerosols by use of the TOMS aerosol index (AI) technique [ Herman 
et al., 1997 Torres et al., 1998]. Specifically, the AI is zero for water clouds, positive for 
absorbing aerosols, and negative for non-absorbing aerosols (e.g., sulfate aerosols). The 
AI is available from standard TOMS data, as a byproduct of the total ozone retrieval 
algorithm. The AI data can be used directly to correct surface irradiance for attenuation 
by absorbing aerosol plumes, if the aerosol altitude is known from either direct 
measurement or climatological estimates: 

17 

1 aerosol _ „ - g (H A ) AI 

F ~ e (l) 

where the conversion factor g is a function of aerosol height, Ha [. Krotkov et al., 1998; 
Herman et al, 1999]. One method for obtaining Ha is by using the GSFC wind-data 
assimilation model to estimate aerosol plume heights by matching modeled and observed 
plume trajectories over several days. This method, and comparisons with sunphotometer 
data, show that the aerosol plumes are usually located between 3 and 4 km altitude, and 
are repetitive at a given location from year to year. Error estimates for uncertainties in 
Ha are given in Herman et al. [1999]. 


3. Model 

The atmospheric RT model provides the boundary conditions at the ocean surface for the 
underwater irradiance calculation. The radiative transfer solutions in the atmosphere and 
in the ocean are coupled through the contribution of photons first reflected from the 
ocean and then scattered back to the water by the atmosphere. The ocean contribution to 
the downward clear-sky irradiance, F c i ear (A), can be estimated from the following 
equation: 


Fciear a) = F 0 (A) (1+T)/ (1 -A S b (A) ) (2) 

where A is the ocean albedo, Fo is direct irradiance at the surface for A=0, r is the ratio 
of diffuse to direct irradiance at the surface for A=0 and Sb(A) is the fraction of reflected 
radiation backscattered to the ocean by the atmosphere [Herman et al., 1999]. /"and S b (A) 
were calculated for a clear atmosphere that contains Rayleigh scatterers and ozone using 
DISORT code [Stamnes et al., 1988]. Given the maximum value of S),=0.4 for a clear 
sky and ocean albedo values less than 10% in the near UV region [Herman and Celarier, 
1997], one can neglect the coupling of the atmospheric and oceanic radiative transfer as a 
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first approximation. We estimate that the flat ocean surface model and the separation of 
the atmospheric and oceanic RT models both give less than 10% resulting error for 
satellite estimation of the underwater UV irradiance. 

3.1 Surface UV calculation 

The combined TOMS datasets, along with the satellite measured high-resolution 
extraterrestrial solar irradiance spectra, are used as an input to our atmospheric radiative 
transfer model to calculate surface UV irradiance at the time of the satellite overpass 
(approximately 1100 LT for Earth-Probe / TOMS). Using daily measurements of ozone 
from TOMS, equation (2) is applied to calculate F c i ear at the surface for clear-sky and 
aerosol-free conditions. The direct radiation is calculated from Beer’s law with airmass 
correction for high solar zenith angle (0 < 75°). The parameters r and Sb are interpolated 
from a look-up table pre-calculated for a Rayleigh atmosphere using climatological 
TOMS 325DU ozone and temperature profiles for different solar zenith angles, assuming 
an average ocean albedo of 6% at 360nm. The angular distribution of the downward 
radiance at the air-ocean boundary is assumed isotropic with a flat ocean surface (Fresnel 
refraction). The absolute irradiances were obtained using extraterrestrial solar irradiance 
data, acquired by the SOLSTICE [Rottman et al., 1993] instrument on NASA's Upper 
Atmospheric research Satellite (UARS). 

A simple model to account for reduction of F c i ear caused by clouds over a low reflecting 
ocean surface is based on an energy balance between solar radiation reaching the ground 
and backscatered to space. For this model [ Eck et al., 1995; Herman et al., 1996; Herman 
et al., 1999]: 


F cloud ~ F clear Cj 

C T ~ 1 -(R 3 60-0.06)/0.88 for 0.06<R 360 <0.5 (3) 

Ct—1 - R 3 60 for R 3 6o >0.5 


This equation accounts for UV irradiance reduction from both clouds and non-absorbing 
aerosols, but needs an additional correction in the presence of absorbing aerosol plumes, 
where UV irradiance reduction is stronger. When plumes of absorbing aerosols are 
detected in the TOMS Al data, the UV irradiance is corrected using Equation (1) 
[Krotkov et al., 1998]. Typical levels of UV reduction can be up to 50% under Saharan 
dust and the annual smoke plumes from Africa and South America. This reduction can 
reach 90% for dense smoke near the fire source and for unusual conditions such as 
occurred in Indonesia and Mexico during 1998 and early 1999 [Herman et al., 1999]. 

The daily UV dose is calculated by integrating over time and assuming that Ct is constant 
throughout the day. This daily cloud variability error is assumed random, therefore it is 
reduced with weekly averaging of daily exposures. Figure 5 shows the 10-day average 
TOMS daily erythemal [McKinlay and Fiffey, 1987] exposure map including the 
combined effects of ozone, clouds and absorbing aerosols [Herman et al., 1999]. The 
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most pronounced effect is the dependence on solar zenith angle, 0 so that geographical 
areas whose latitudes are near sub-solar points (e.g., +23N for July) produce very high 
clear-sky surface UV irradiance. The next strongest source of surface UV modulation is 
from clouds. In July, there is considerable haze and clouds over the ocean areas in the 
Northern Hemisphere, since the Intertropical Convergence Zone (ITCZ) cloud region has 
moved toward the north. In terhemi spheric differences in cloud cover cause the Southern 
Hemisphere summer to have even higher exposures than the Northern Hemisphere 
summer for the same SZA [Herman et al., 1999]. The remaining geographical variations 
can be attributed to ozone variations and absorbing aerosols in certain regions. 

To provide boundary condition to the underwater RT, we must separate the TOMS 
estimated global irradiance (F c i 0M j) into direct and diffuse components. The direct 
component is estimated using the SeaWiFS cloud fraction data, Scf (in percent): 

F direct = Fo (A) ( 1 - 0.01 Scf) (4) 

Finally, the diffuse component at the surface is determined as the residue between the 
TOMS global irradiance, Fci 0U d, and the direct component estimated from SeaWiFS cloud 
fraction in equation 3: 

F diffuse — F cloud ~ F direct (5) 


3.2. Radiative transfer in the ocean 

Given the TOMS estimate of the surface UV irradiance, and assuming isotropic angular 
distribution of the diffuse downward radiance at the ocean surface, any appropriate 
radiative transfer scheme can be applied to model light penetration into the ocean [Jin 
and Stamnes, 1994; Mobley et al., 1993]. However, the accuracy of the current optical 
measurements of the most fundamental inherent optical properties of seawater (scattering 
and absorption coefficients) is normally about 10%, and the errors of extrapolation of 
these properties into UV spectral region has yet to be estimated. This currently prevents 
us from applying these sophisticated radiative transfer schemes for the purpose of 
operational satellite mapping of underwater UV fields. However, the accurate models 
[Mobley et al., 1993] are very important for testing faster operational algorithms. 

3.2.1. The Quasi-single scattering approximation of radiation transfer in water 

The first model for an assessment of underwater UV radiation and biodoses was 
developed by Smith and Baker [1979]. In this model it was assumed that the irradiance is 
attenuated exponentially in the ocean: E\(z) = E° exp(- K d z), where Ex and Ex(z) are 
spectral irradiances just below the ocean surface and at depth, z (we use different notation 
E for irradiances within the ocean to distinguish them from the atmospheric irradiances: 
F), and the diffuse attenuation coefficient, K d , is assumed to be vertically constant. A 
similar model was applied by Booth et al. [1992]. This simple formulation of the 
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radiative transfer in the ocean requires an a-priori knowledge of Kd in the UV spectral 
region. 

In general, Kd cannot be extrapolated from the visible region for use in the UV 
wavelengths. The coefficient Kd depends on the angular structure of the light field and, 
thus, on depth (even for a homogeneous ocean), and on seawater inherent optical 
properties (IOPs). Therefore, there is no a-priori reason to expect that Kd values in the 
UV region will vary in the same manner with the angular structure of the light field and 
depth as in the visible region. However, in some ocean areas this correlation was 
reported. For example, Hojerslev and Aas (1991) have shown that ^(310) is linearly 
correlated to Kd(465) for vast ocean areas with the correlation coefficient as large as 
0.998. The problem of correlation between spectral values of the diffuse attenuation 
coefficient is carefully discussed by Hojerslev [1982, 1986]. From the theoretical point of 
view the relationship between the diffuse attenuation coefficient and seawater IOPs is 
considered by Gordon [1989]. 

To calculate UV underwater irradiances, we chose here the, so-called, Quasi-Single 
Scattering Approximation (QSSA) [ Gordon , 1973]. The QSSA model is important, 
because of its simple analytical formulation, yet enabling us to address the dependence of 
Kd on the angular distribution of the light field in the ocean. Thus, the model can account 
for Kd dependence on depth even in vertically homogeneous water. The QSSA is based 
on the strong absorption with highly anisotropic scattering of seawater [Gordon, 1973]. It 
assumes: (a) single scattering in the upward direction; (b) multiple scattering in the 
downward direction in accordance with delta-function. The approximation assumes the 
exponential direct-beam transmittance: T b = exp( -Kz/ju), where K = a + b\ is the 
effective attenuation coefficient, a is the absorption coefficient, b \ is the fraction of 
backscatter in the upward direction, and ju is the cosine of the zenith angle of incident 
beam. The fraction of backscatter in the upward direction is exactly equal to the 
backscattering coefficient, bb, in the case of normal incident beam, i.e. ju-l. In-water 
zenith angles are limited by the angle of total internal reflection. Accounting for this fact, 
and the small contribution of backscatter to total attenuation of light, we assume that 
b b=h. 

Here the water column is assumed to be vertically homogeneous. This assumption can be 
easily accepted because the penetration depth for UV radiation does not ordinary exceed 
the depth of the upper quasi-homogeneous, well-mixed layer in the ocean. The spectral 
irradiance at depth, z, can be written as a sum of the direct solar radiation and the integral 
of surface radiance over spherical angles from the diffuse radiation, both attenuated by 
water [Krotkov and Vasilkov, 1993]: 


2n Ml 

Ex(z) = Ex T b (Mo) + /d<p/U(M,<P) T b (n) M dM (6) 

0 0 

where E\ is the spectral irradiance for direct solar radiation just beneath the ocean 
surface, Mo is the cosine of the solar zenith angle in the water, (p is the azimuth angle, M\ is 
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the cosine of the total internal reflection angle -49°, L\(p,(p) is the angular distribution 
of sky radiance after surface refraction and reflection. 

3.2.2. Model of seawater inherent optical properties (IOPs) 

The QSSA makes use of the absorption coefficient, a , and the backscattering coefficient, 
b b . The total IOPs are the sums of the IOP of the pure seawater and the three major 
scattering and absorbing water substances: 

a(A) = a w (A) + odom(A) + ciph(A), bt(A)= bw(A) + b p (A) (7) 

where subscripts w, p, ph and DOM denote the pure seawater, the particulate matter, the 
phytoplankton pigments, and dissolved organic mater (DOM), respectively. The detritus 
absorption is included into the DOM absorption because of approximately identical 
spectral dependence [Carder et al., 1991]. 

The pure seawater IOPs were obtained from Smith and Baker [1981]. According to the 
recent findings [Pope and Fry, 1997; Sogandares and Fry, 1997] the pure water 
absorption coefficient is significantly below the consensus values (as previously selected 
by Smith and Baker [1981] ) in the wavelength range 380 - 500 nm, about 2 times lower 
than the old value at 380nm. Currently, no revised measurements of the pure-water 
absorption coefficients are available for wavelengths shorter than 340nm. Therefore, for 
the present paper we use values of the pure-seawater absorption coefficient from Smith 
and Baker [1981] in most calculations. 

The suspended particulate matter (SPM) backscattering coefficient and the DOM 
absorption coefficient are used in the conventional form: 

odoM = a 0 expl -S(A-U)] , b b (A)= bd( W" (8) 

where m is the backscatter wavelength ratio exponent, S is the DOM spectral slope. 

The DOM spectral slope 5=0.014 nm' 1 is commonly accepted for the visible spectral 
region [Bricaud et al., 1981]. A study by Kopelevich et al. [1989] showed that the DOM 
spectral slope should be put slightly greater in the 280-490 nm spectral region: 5=0.017 
±0.001 nm' 1 . More recent measurements have showed that the DOM spectral slope can 
increase with photodegradation of colored DOM and can vary within a rather wide range 
from 0.01 to 0.03 nm' 1 for clear waters [Vodacek and Blough, 1997]. Unfortunately, these 
variations of the DOM spectral slope have not been parameterized. We use, therefore, an 
average value of the DOM spectral slope for the UV spectral region 5=0.017 nm' 1 . The 
parameter m may vary in a wide range depending on the optical type of seawater. 
Fortunately, the SPM backscattering coefficient, bo, is normally much less than the total 
absorption coefficient, a, in the UV spectral region. Therefore, an average estimate of the 
parameter m - 1 was adopted in the model of Gordon and Morel [1983]. 
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The phytoplankton pigment absorption is expressed through chlorophyll a concentration 
and the chlorophyll-specific absorption coefficient: 


a Ph (A) = Ca* Ph (A) (9) 

Phytoplankton synthesize a variety of compounds that absorb radiation in the UV-B and 
UV-A regions of the spectrum and which could affect the response of the cell to UV 
radiation. In vivo absorption in the UV-A and UV-B shows a wide range of values, with 
peaks of absorption between 320 and 350 nm and a maximum between 330 and 335 nm 
(Vemet et al., 1994). For all surface samples, the average in vivo pigment-specific UV 
absorption was larger than absorption in the blue. The overall variance was high, 
demonstrating that the absorption in the UV is not due to the major photosynthetic 
pigments and that the UV-absorbing compounds, such as mycosporine-like amino acids, 
vary independently. The average for all-stations spectrum was accepted to be the 
chlorophyll-specific absorption coefficient in the model. Values of the function 
R(A)=a P h(A) / a P h (440) calculated on basis of this data are given in Table 1. Values of 
the function calculated in [Kopelevich and Ershova, 1997] are also given in Table 1 for 
comparison. As can be seen from Table 1 the values A(A) [Kopelevich and Ershova, 
1997] are less than the values R(A) [Vemet et al., 1994] even for the average minus 
standard deviation spectrum. We use the data of Vemet et al. [1994] in our calculations. 

Table 1. The function R{A) calculated on basis of data of [Vemet et al., 1994] 


Wavelength, nm 

380 

360 

340 

331 

Average plus standard deviation spectrum 

0.81 

0.90 

1.56 

1.80 

Average spectrum 

0.78 

0.80 

1.27 

1.47 

Average minus standard deviation spectrum 

0.65 

0.60 

0.65 

0.85 

after Kopelevich and Ershova, 1997 

0.49 

0.42 

0.61 

0.68 


The model contains three input quantities: ao, bo, and C. These parameters are to be 
estimated from available satellite data sets. Some seawater characteristics can be 
retrieved in the visible spectral region using the data of satellite ocean color sensors. The 
problem of first importance is an extrapolation of the water absorption and scattering 
coefficients measured or retrieved from satellite measurements in the visible (400 - 600 
nm) into the UV spectral region (290 - 400 nm). The above spectral model of IOPs 
allows such an extrapolation provided the model-input parameters are known and 
parameterization of IOPs is valid. 

What should be definitely refined in the future is the parameterization of phytoplankton 
absorption. It is well known that the chlorophyll-specific absorption coefficient depends 
on chlorophyll concentration due to, for example, pigment packaging effect. This 
dependence has been parameterized for the visible range in Bricaud et al. [1995]. 
However, for the UV region, the more important issue may be a weak correlation 
between photosynthetic pigments and UV-absorbing, photoprotective pigments. Because 
of lack of such data, we chose a rather simple parameterization of the chlorophyll- 
specific absorption coefficient. 
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Next, we will consider how the model-input parameters can be estimated from available 
standard products of satellite ocean-color sensors. First, the chlorophyll concentration is 
the standard SeaWiFS product. To determine other quantities, the Case 1 water model 
[Morel, 1988] is assumed. According to the model, the DOM absorption at 440 nm is 
20% of the total absorption of pure seawater and pigments. This assumption determines 
the most important parameter ao. To estimate the backscattering coefficient we use the 
standard SeaWiFS product of the diffuse attenuation coefficient [McClain et al., 1992] 
and the model of the diffuse attenuation coefficient developed in Gordon [1989]: 

K d /D 0 = 1.0395 (a + b b ), (10) 

Do =f/cos0o„ +1.197 (l -f), 

where /is the fraction of direct sunlight in the incident irradiance transmitted through the 
sea surface. The fraction /is calculated as a function of solar zenith angle and wavelength 
from the atmospheric model. 

Coastal waters are normally referred as Case 2 waters in which IOPs are uncorrelated. 
The above model cannot be directly applied to those waters. However, independent 
retrieval of absorption coefficients of DOM and phytoplankton pigments has been 
suggested in Carder et al, [1991]; Hoge and Lyon, [1996]; Garver and Siegel, [1997] 
and Vasilkov, [1997]. Given the DOM absorption coefficient in the visible region, it can 
be extrapolated into the UV region. The SPM backscattering coefficient can also be 
retrieved using analytical algorithms [Hoge and Lyon, 1996', Garver and Siegel, 1997; 
Vasilkov, 1997]. 

3.2.3. Verification of the QSSA 

The accuracy of the QSSA was estimated by comparison with the accurate RT 
calculations for the simplified models of the ocean. In this paper we used radiative 
transfer code intercomparisons [Mobley et al., 1993] and results of an accurate Monte- 
Carlo radiative transfer scheme [Vasilkov et al., 1999]. The underwater irradiance was 
compared with data for Problem 2 from the paper by Mobley et al., [1993]. The same 
scattering phase function was used in the QSSA calculation. The relative error of the 
QSSA within the optical depth range of 10 was less then 7% for the single scattering 
albedo &*=0.2. For the single scattering albedo ft*=0.9, the QSSA error was 17% for the 
optical depth and 49% for the optical depth £=10. 

Problem 2 [Mobley et al., 1993] deals only with the direct solar radiation illuminating the 
surface. To estimate the errors of the QSSA in case of the diffuse illumination, 
Monte-Carlo calculations were conducted for an isotopic angular distribution of incident 
radiance on the ocean surface [Vasilkov et al., 1999]. For a typical single scattering 
albedo to=0.6, the QSSA error was less than 35% in the optical depth region t< 10. In all 
cases the errors were less at smaller optical depths. Similar results regarding the QSSA 
accuracy has been reported using Monte-Carlo simulations [Greysukh and Levin, 1990]. 
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It should be noted that certain conditions specific for UV spectral region justify the 
application of the QSSA to the radiative transfer problem in the water: 

(a) to is normally less than 0.7 in the UV spectral region and gets smaller at short 
wavelengths, which are more biologically effective; 

(b) only small optical depths play a significant role in biological applications of the 
underwater UV calculations. All these considerations should permit the use of the 
QSSA for calculations of biologically significant parameters from the underwater UV 
irradiance. 

4. Results 

The biological effect of UV radiation is typically described by action spectra. A large 
number of action spectra, A(X), has been proposed for various biological effects of UV 
radiation in marine environment (for Antarctic waters see discussion by Smith et al. 
[1992]). Two action spectra with considerably different shapes have been chosen for a 
sensitivity study: (1) the marine phytoplankton [Mitchell, 1990] and (2) unshielded DNA 
(Setlow spectrum) [Quaite et al., 1992]. We simulate the typically measured biological 
daily UV doses by convolution of UV irradiance spectra- Exiz) with A(A) and integrating 
over the time of the day: 

400 

D(z) =fdt fExizMt)) A(A) dA (1 1) 


To estimate the damage for particular biological process, the units of Ex should agree 
with the units of A(A). For example, to model the conditions of the field experiments 
described in Reagan et al. [1992], we have chosen the action spectrum for unshielded 
DNA [Quaite et al., 1992]. Since the A(A) is given in terms of number of pyrimidine 
dimers per photon/m , the value of Ex in W / (m nm) was transformed to quantum units 
as E p h -Ex/ he, where h is the Plank constant and c is the light speed. The calculations of 
D(z=0) gave a value of 0.5% pyrimidine dimers for typical total column ozone values of 
260DU near the Bahamas (23°45’N) in July 1990. This is close to the estimate by Reagan 
et al. [1992], that daily exposure of DNA molecules just above the surface produced 
0.45% pyrimidine dimers (or 4.5 dimers per thousand bases of molecule). Unfortunately, 
no water optical properties were measured in the field experiment by Reagan et al. 
[1992], which prevents a comparison with underwater data. 

4.1 DNA doses 

The weekly averaged DNA dose at depth of 3 m is shown in Fig. 6. It is seen that the 
main features of the averaged DNA dose map are determined by latitude dependence of 
the surface UV irradiance. The latitude dependence of the DNA dose is clearly apparent 
in all oceans. Some features of the DNA dose map are due to cloudiness structure. For 
example, reduced levels of the DNA dose in some areas of the Indian and Pacific Oceans 
correspond to 100% cloudy skies observed in these areas (see map of cloud fraction in 
Fig. 4). The cloudiness effect on the DNA dose is observed in the Mediterranean Sea - 
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where clear-sky conditions remained for more than a week, resulting in DNA dose values 
characteristic of equatorial regions. 

It is interesting that the latitudinal distribution of both total ozone and the optical 
properties of ocean waters are not seen on the global DNA dose map. The effects of 
ozone amount and seawater optical properties are almost masked by cloudiness effects, 
indicating that latitudinal dependence of the UV irradiance and cloudiness are the major 
factors affecting the underwater DNA dose. 

4.2 Z10 map 

An important measure of the DNA dose is the depth at which DNA dose is reduced to 
10% of its surface value (or Z10). This is the approximate depth over which biological 
damage due to UV effects take place. The 10% depth depends somewhat on the action 
spectrum used in calculations of UV dose rates. Since this dependence is rather weak, the 
10% DNA dose depth can be used to characterize the biologically significant depth range. 

A global map of the 10% DNA-dose depth is shown in Fig. 7a. Comparing this map with 
the chlorophyll and A^(490) maps, one can easily see that horizontal distribution of the 
10% DNA dose depth is primarily determined by bio-optical properties of ocean waters. 
However, the angular structure of the light incident on the sea surface also affects the 
10% DNA dose depth. This is because the dependence of the diffuse-attenuation 
coefficient on the angular structure of the in-water light field. The angular structure of the 
light incident on the sea surface is determined by the SZA and cloudiness. Therefore, a 
slight latitude dependence and cloudiness structure can be seen in the map of the 10% 
DNA dose depth. 

A global map of Z10 is shown in Fig. 7b for unweighted UVB irradiance (280nm- 
320nm). The 10% UVB irradiance depth is normally greater than Z10 (the 10% DNA- 
dose depth). This is due to the fact that the integral over UVB irradiance is mainly 
determined by the longer wavelength part of UVB spectrum as opposed to the DNA dose 
that is mainly determined by shorter UVB wavelengths. The seawater is an effective filter 
of the shorter UV wavelengths. 

Z10 is mainly determined by seawater IOPs rather then by the angular distribution of 
light incident on the sea surface. We carried out a sensitivity study intended to estimate 
the errors in Z10 due to uncertainties of the seawater IOPs. 

The absorption coefficient of pure seawater has a strong affect on the UV penetration 
depth. To estimate the effect of its uncertainty on the UV penetration depth, data on the 
water absorption coefficient available in the spectral region of 340-380 nm [Sogandares 
and Fry, 1997] was extrapolated to 290-340 nm by using a least-squares technique 
applied to logarithms of the water absorption coefficient. Linear extrapolation resulted in 
lower values of the water absorption coefficient than logarithmic extrapolation. No 
absorption by salt, Chi a, or DOM was assumed. A comparison of the diffuse attenuation 
coefficient and penetration depths calculated for different water absorption coefficients is 
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shown in Figure 8, along with the diffuse attenuation coefficient measured in clear ocean 
waters [Booth et al., 1992]. 

The Z10 UVB penetration depth calculated from the water absorption coefficient of 
Sogandares and Fry [1997] is about 20% greater than that calculated from coefficients of 
Smith and Baker [1981]. For example, the penetration depth at 320 nm is 29.9 m using 
the Sogandares and Fry [1997] data and 23.3 m for the Smith and Baker [1981] data. The 
10% UVB depths are 24.2 m and 19.9 m respectively. It should be noted that the diffuse 
attenuation coefficients measured at 308, 320, and 340 nm are only slightly greater than 
those calculated using the Sogandares and Fry [1997] data. Booth and Morrow [1997] 
compared various measurements of the UV penetration in natural water and found no 
penetration depths at 320 nm greater than about 23m. This means that the 10% UVB 
irradiance depth generally should be less than 23 m. The only exception is data recently 
presented by Aas and Hojerslev [1999]. They reported of the 10% UVB irradiance depth 
as large as 35 m. The sensitivity study highlights the importance of accurate knowledge 
of the pure water absorption coefficient. 

It is instructive to estimate how variations in the chlorophyll concentration and diffuse 
attenuation coefficient affect the UVB penetration depth. A contour plot of depths of 
UVB irradiance penetration to a level of 10% of its surface value is shown in Figure 9. In 
this case, the UVB penetration depth was calculated for the seawater absorption 
coefficient after Smith and Baker [1981]. No values of the UVB penetration depth for 
low 7fd(490) are available because the necessary condition Kj > a + bb is not valid 
[Gordon, 1989]. Calculations were conducted for two cases. In the first one, no CDOM 
absorption was assumed. The case represents upper limit values of the penetration depth. 
In the second case, it was assumed that the DOM absorption at 440 nm is 20% of the total 
absorption of pure seawater and pigments [Morel, 1988]. Figure 9 demonstrates the 
significant effect of the DOM absorption on the UVB penetration depth. For example, the 
UVB penetration depth for most probable chlorophyll concentration of 0.1 mg/m 3 and 
Kd{ 490) = 0.03 m' 1 decreases from 17.4 m to 14.1 m. The DOM absorption reduces the 
maximum penetration depth from 19. 9m to 15. 8m. 

5. Conclusions 

A radiative transfer model for the atmosphere-ocean system has been developed for 
purpose of mapping underwater UV irradiance and biologically active UV dose rates in 
the ocean. The model is intended to combine total ozone and UV reflectivity 
measurements from TOMS, with ocean-color data available from CZCS and SeaWiFS. 
The ocean-color data includes monthly maps of chlorophyll concentration and the diffuse 
attenuation coefficient as well as daily maps of cloud fraction averaged on a grid of 1 by 
1 degree. 

The atmospheric part of the model generates spectral direct and diffuse irradiance on the 
sea surface that are inputs to the underwater part of the radiative transfer model. The in- 
water radiative transfer model is based on the QSSA (Quasi Single Scattering 
Approximation) that is simple, computationally fast, and yet enables the angular 
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distribution of the light field to be addressed. A model of seawater IOPs allows the 
extrapolation of the absorption and backscattering coefficients to the UV spectral region 
provided their values in the visible region are known. Values of the absorption and 
backscattering coefficients in the visible region are estimated from the SeaWiFS standard 
products by using the Case 1 water model. 

The sensitivity study has shown that the main parameters controlling levels of the most 
harmful UV-B radiation underwater for clear sky conditions are the solar zenith angle, 
water bio-optical properties and total ozone. Attenuation of UV-B irradiance and DNA 
dose rate with water depth is primarily controlled by the seawater absorption coefficient 
and its spectral dependence. An influence of the seawater backscatter on the attenuation 
of UV irradiance is considerably less. Changes in the angular distribution of the surface 
radiance due to aerosol load or clouds may result in the irradiance increase at a given 
depth for large solar zenith angles. The sensitivity of UV-B and DNA dose rate to ozone 
variations decreases with depth. 

The main spatial features of the weekly maps of underwater DNA dose are determined by 
the SZA and cloudiness. The seawater IOPs and total ozone effects are less significant for 
the spatial distribution of the DNA dose. The spatial distribution of the 10% DNA dose 
depth is mainly determined by the spatial structures of chlorophyll. Cloudiness effects 
and latitude dependence of the 10% DNA dose are also observed due to the effect of the 
angular distribution of the light incident on the sea surface in the in-water UV irradiance 
attenuation 
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Figure Captions 


Fig. 1. Data flow chart showing how TOMS and SeaWiFS products are assimilated by 
the suggested model. 

Fig. 2. SeaWiFS-derived monthly chlorophyll a map for July 1998. 

Fig. 3. SeaWiFS-derived monthly ^(490) map for July 1998. 

Fig. 4. Map of the SeaWiFS-derived cloud fraction averaged over 5 days. 

Fig. 5. TOMS-derived map of 5 day averaged surface UV irradiance from 290 to 400 nm 
weighted by the erythemal action spectra [McKinlay and Dijfey, 1987]. The energy is 
expressed in kilojoules per square meter. 

Fig. 6. Weekly averaged DNA dose map for July 1998. The DNA dose is expressed in a 
number of pyrimidine dimers produced per thousand bases of the DNA molecule during 
daytime. 

Fig.7a. Map of the weekly averaged 10% DNA dose depth for July 1998. 

Fig.7b. Map of the weekly averaged 10% UVB flux depth for July 1998. 

Fig. 8. The diffuse attenuation coefficient and depth of penetration to level of 10% of the 
surface irradiance. Solid lines represent results calculated for the water absorption 
coefficient after Smith and Baker [1981]. Dashed lines represent results calculated for the 
extrapolated water absorption coefficient after Sogandares and Fry [1997]. Asterisks 
denote values measured for clear ocean waters [Booth et al., 1992]. 

Fig. 9. The depth of UVB penetration to level of 10% of the surface irradiance as a 
function of the input parameters: chlorophyll concentration and the diffuse attenuation 
coefficient at 490 nm. (a) No DOM assumed, (b) DOM assumed to be according with the 
Case 1 water model [Morel, 1988]. 
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Popular Abstract 


GLOBAL MAPPING OF UNDERWATER UV IRRADIANCES AND DNA- 
WEIGHTED EXPOSURES USING TOMS AND SEAWIFS DATA PRODUCTS 

Alexander Vasilkov 1 , Nickolay Krotkov 1 , Jay Herman 2 , Charles McClain 2 , Kevin 
Arrigo 3 , and Wayne Robinson 4 


Abstract. Global or local increases in penetration of UV radiation in the ocean can have 
harmful effects on the food chain. The global stratospheric ozone-layer depletion results 
in an increase in biologically harmful ultraviolet (UV) penetrating to ecologically 
significant depths in the ocean. Such an increase can be estimated anywhere in the 
Earth’s oceans by combining satellite estimates of UV surface irradiance from the Total 
Ozone Mapping Spectrometer (TOMS) satellite instrument with the SeaWiFS satellite 
ocean-color measurements. We use a model of seawater optical properties with input 
from standard monthly SeaWiFS products: chlorophyll concentration and the diffuse 
attenuation coefficient at 490nm. Penetration of solar UV radiation to different depths in 
open ocean waters is calculated using an approximation to the full RT (radiative transfer) 
calculation with the results compared to more accurate models. We show that the the 
main environmental parameters controlling the absolute levels of the UVB (280-320nm) 
and DNA-weighted irradiance underwater are: solar-zenith angle, cloud transmittance, 
water optical properties, and total ozone. Weekly maps of underwater UV irradiance and 
DNA-weighted exposure are calculated using monthly-mean SeaWiFS chlorophyll and 
diffuse attenuation coefficient products, daily SeaWiFS cloud fraction data, and the 
TOMS-derived surface UV irradiance daily maps. The final products include global maps 
of weekly- average UVB irradiance and DNA-weighted daily exposures at ocean depths 
of 3m and 10m, and the depths where the UVB irradiance and DNA-weighted dose rate 
at local noon are equal to 10% of their surface values. 
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Monthly Average SeaWiFS Chlorophyll Concentration, July 1 990 


Figure 2 


Monthly Average SeaWiFS Kd(490nm), July 1 99B 
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Figure 3 



Average SEAWIFS Cloud Fraction July 14 through 19, 1998 
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DNA Exposure at 3m [Dimers/1000] July 9 to 19, 1998 


Figure 6 



Average Z10Z DNA [meters] July 9 to 19, 1998 


Figure 7a 




Average Z10% UVB [meters] July 9 to 19, 1998 


Figure 7b 











